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Abstract. Mode conversion in the region where the sound and Alfven speeds are equal is a 
complex process, which has been studied both analytically and numerically, and has been seen 
OO ■ in observations. In order to further the understanding of this process we set up a simple, one- 
dimensional model, and examine wave propagation through this system using a combination of 
analytical and numerical techniques. Simulations are carried out in a gravitationally stratified 
atmosphere with a uniform, vertical magnetic field for both isothermal and non-isothermal cases. 

(3JT)' For the non-isothermal case a temperature profile is chosen to mimic the steep temperature 
gradient encountered at the transition region. In all simulations, a slow wave is driven on the 
■ upper boundary, thus propagating down from low-/? to high-/? plasma across the mode-conversion 
region. In addition, a detailed analytical study is carried out where we predict the amplitude and 
phase of the transmitted and converted components of the incident wave as it passes through 
the mode-conversion region. A comparison of these analytical predictions with the numerical 
results shows good agreement, giving us confidence in both techniques. This knowledge may be 

O ,■ used to help determine wave types observed and give insight into which modes may be involved 

^ ' in coronal heating. 

' Keywords. Sun: oscillations, MHD, gravitational waves. 
,— 1 ' 1. Introduction 

ly-j ■ Mode conversion is the process under which an incident wave may be wholly or par- 
, tially converted into another wave mode, with the remainder of the original wave being 
transmitted through the plasma. This may occur where the sound and Alfven speeds are 
equal, or equivalently, where the plasma /? (the ratio of the gas to the magnetic pressure) 
is approximately unity. This layer will generally lie low down in the solar atmosphere, in 
, the region of the lower chromosphere. 

This problem has been investigated extensively in the past, many using analytical tech- 



niques. Zhugzhda (1979) Zhugzhda fc Dzhalilov (1981) and Zhugzhda fc Dzhalilov (1982) 



investigated the conversion of fast magnetoacoustic waves into slow magnetoacoustic 
waves as they propagate up through the solar atmosphere from high- to low-/? plasma. 
Their model looked at a uniform vertical magnetic field set in an isothermal atmosphere 
stratified by gravity. Using the fact that, for harmonic waves, an exact solution may be 
found in terms of the known Meijer-G functions, coefficients describing the conversion 



and transmission were found. Cally (2001) expanded on this work by noting that these 



solutions may also be written in terms of the equivalent Hypergeometric 2^3 functions, 
and another set of coefficients was discovered. 

We also investigate the mode conversion problem analytically, although we combine 
this with numerical simulations. We look at a uniform background magnetic field, but, 
in contrast to the aforementioned work, we investigate downward propagation - so the 
wave is travelling from low- to high-/? plasma (as shown in Figure [J) . As the fast wave 
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Figure 1. Cartoon of our model atmosphere, with uniform vertical magnetic field. The z-axis 
points upwards (opposite to gravity) and a slow wave is driven on the upper boundary, travelling 
down towards the mode-conversion layer at /3 « 1. 

is evanescent in low-/? plasma, we look at the conversion of slow to fast magnetoacoustic 
waves. 

In Section [2] the basic equations are set out, and we look at the characteristic wave 
behaviour in high- and low-/? plasma. Two different situations are then examined: an 
isothermal atmosphere in Section [31 and a non-isothermal atmosphere in Section [U Fi- 
nally the conclusions are detailed in Section [5] 



2. Basic Equations 

Throughout we use the standard, ideal MHD equations, linearised about the equilib- 
rium Bo = -Bok, and p' = —pog', where po = PqRTq (z) and H = g/RTo may be defined 
as the scale height. So we have the equation of motion 

Po-J 1 = -Vpi + -(VxB 1 )xB + p 1 g, (2.1) 

Ot p 



the induction equation 

the mass continuity equation 

and the energy equation 



ggi 
~df 

dpi 
dt 



V x (vi x B ) , (2.2) 
V-(p o vi) = 0, (2.3) 



where subscript zero denotes the equilibrium quantities and subscript one the perturbed 
quantities. In these equations B is the equilibrium magnetic field, and po and p are the 
equilibrium pressure and density which may be calculated from the equilibrium equation 
above. The variable vi denotes the plasma velocity, p\ is the plasma pressure, p is the 
magnetic permeability, Bi is the magnetic field, p\ is the mass density, g = —5k is the 
gravitational acceleration, and 7 is the ratio of specific heats, which we take to be 5/3. 
Henceforth, the subscripts on the perturbed variables are dropped, and it is understood 
that we are working with the linearised equations. 

Under the assumption that the x-dependence is of the form expik x x, where k x is 
the horizontal wavenumber, Equations (|2.ip - (|2 .4[) may be combined to form a pair of 



298 



A. M. D. McDougall and A. W. Hood 



wave equations 



i (cl + v\) k 2 x - ^\ v x = k x (c 2 s ^- -g)v z , (2.. r > i 



)z 



d 2 d d 2 \ ( d \ 

cl dz- 2 ~ 19 d- z ~dT 2 ) Vz = ~ kx [ cl d- z ~ (7 - 1)9 J (2 ' 6) 

where c 2 — jpo/po and v\ = B 2 / (fipo) are the squared sound and Alfven speeds respec- 



tively. These are in agreement with those discovered by Ferraro & Plumpton (1958) 

Examining these equations we can see that if the effects of gravity are negligible and 
the horizontal wavenumber is small, then the wave operators on the left-hand side are the 
same when c s = va- At this point there is a resonance which allows the mode conversion 
process to occur. Note that when the horizontal wavenumber is zero, there is no mode 
conversion. 

We solve these equations numerically using the MacCormack finite difference scheme. 
This is a type of Lax-Wendroff, two-step, predictor-corrector method. We use backward 
differencing for the predictor steps and forward differencing for the corrector steps, thus 
ensuring that we have more accurate values on the upper boundary where we drive a 
sine wave in the vertical velocity v z . As previously mentioned, this is predominantly a 
slow wave as the fast wave is evanescent in the low-/? plasma. 
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Fast Wave 


Speed~ c s 


Speed~ va 




Approx. isotropic 
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Table 1. Characteristic properties of the slow and fast magnetoacoustic waves depending on 

their position in the solar atmosphere. 

As the wave propagates across the /? ~ 1 layer it is important to remember that 
the slow and fast waves have different properties above and below this region; Table [T] 
identifies these. We are driving a slow wave in low-/? plasma, the transmitted component 
of this wave will demonstrate the same characteristics in the high-/? plasma. From the 
table we see that the transmitted wave is thus the fast wave, and the converted part 
must then be the slow wave. 



3. Isothermal Atmosphere 

To begin with we investigate an isothermal atmosphere, so the sound speed will remain 
constant while the Alfven speed varies with height. 

3.1. Numerical Simulations 

The left-hand plot in Figure [5] shows the vertical velocity from the numerical simulation. 
This displays a strong exponential nature which is due to the gravitational stratification, 
and this masks what is occurring at the mode-conversion region (dashed line) which lies 
at z = 0. We may overcome this by making the transformation v z — > v z e z l 2 so that 
the exponential behaviour is removed and the incident wave has a constant amplitude. 
Upon doing this we obtain the right-hand plot of Figure O There is now a clear change 
at the conversion region, with a definite decrease in amplitude and the presence of two 
wave modes. The transmitted fast wave can be seen propagating out in front, while 
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Figure 2. Left: Vertical velocity, v z , for lo — 2-k\/~& and k x — ir at t = 13.5 Alfven times. 
Right: Transformed vertical velocity, v z , for the same parameters. In both cases the dashed line 
indicates where c s — va- 

the converted slow wave is visible in the form of interference. As both the incident and 
transmitted waves have a constant amplitude after this transformation we may calculate 
the proportion of transmission directly. We shall utilise this later to test the accuracy of 
the analytical predictions. 

3.2. Analytical Study 

The first analytical technique we use to examine mode conversion is the WKB method. 
This involves expanding the horizontal and vertical velocity components, v x and v z , in 
inverse powers of lo (the driving frequency). This is valid under the assumption that lo 
is large, and will give solutions valid outside the mode-conversion region. Substituting 
these expansions into the Wave Equations (|2.5| and (|2.6|) we obtain the solutions 



Low/3: Incident Wave v z cx e z / 2 exp ( ) 

\Cs J 

High/3: Transmitted Wave v z oc Be z ^ 2 exp 



ILOZ 
Cs 



n * jw a ik x c 2 8 v A e z/A ( 2iu 
Converted Wave v z oc —A — — = ^— exp 

where A and B are the conversion and transmission coefficients respectively. Note that 
we have also extracted the exponential dependence seen in the simulations. 

These coefficients are the same as those found previously by Zhugzhda & Dzhalilov, 
which were possible to determine because the exact analytical solution was known. We use 



a method developed by Cairns & Lashmore-Davies (1983) to find the coefficients. This 



method is valid at the mode-conversion region, and has the advantage that coefficients 
may be found even when there is no exact solution. We shall take advantage of this when 
we look at the non- isothermal atmosphere in Section 0J 



From Cairns & Lashmore-Davies (1983) we can write the Wave Equations l]2.5p and (|2.6 



in the general form 



— ^ - i ( k - ~£ ) <j>\ = i\<t>2, (3.1) 
( k - 7O 02 =iX<h, (3.2) 
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Figure 3. Ratio of the slow wave amplitudes in high and low f3, telling us how much of the 
incident wave will be transmitted into the high-/? plasma. Here uj — 4ty\ // 6, the solid line is that 
predicted analytically, and the stars are the values calculated from the numerical data. 

where and 4>2 are the main variables, and the conversion and transmission coefficients 
are given by 



/ 27rA 2 a/\ / 7rA 2 a/ \ , . 

In order to manipulate Equations (|2.5[) and (|2.6[) into the correct format we take Fourier 
components in time, we must then impose the conditions that k x is small and uj is large, 
allowing us to neglect some terms. Finally making the transformation v z — ic s V z /vA and 
expanding z = + £ about the mode-conversion region our equations may be written 

dv x . f uj uj \ ik x 

^-*U-^h = ir^ (3 - 4) 

~Tf ~ —Vz = ^v x . (3.5) 
d£ c s 2 

By comparing these to Equations J3.1[) and (|3.2p we may simply write down the conver- 
sion and transmission coefficients 



A=^-ex P (-^j, B.exp^-^— J. (3.6) 

As the amount of transmission is easily calculated from the numerical simulations we 
may compare this directly to what is predicted analytically by Equation (|3.6[) . We did 
this for various values of k x and the results are shown in Figure [3] where the stars indicate 
the transmission in the numerical simulations and the solid line is that predicted ana- 
lytically. The agreement between the numerical and analytical results is excellent, even 
when k x becomes large thus violating our assumptions. This suggests that the analytical 
predictions are an excellent tool for quantifying mode conversion. More information on 



the analytical techniques may be found in McDougall & Hood (2007) 



4. Non-Isothermal Atmosphere 

Exactly the same analysis may be carried out for a non-isothermal atmosphere, where 
there is no exact analytical solution. We choose a tanh curve for the temperature profile, 
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Figure 4. A temperature profile of the form A = a + fctanhz is chosen to mimic the steep 
gradient of the transition region. Here we have set a = 0.55 and b = 0.45, and the dashed line 
indicates where c s — va- 

as shown in Figure 21 as this mimics the steep temperature gradient of the transition 
region. The form of this profile is A = a + fetanh(z/i?) and c s = va is chosen to lie 
in the centre of the steep gradient at z = 0, as indicated by the dashed line. We may 
use Wave Equations ()2.5() and (|2.6|) here, but must remember that both the sound and 
Alfven speeds will now vary with height. 

4.1. Analytical Study 

Using a WKB analysis on Equations (|2.5p and (|2.6p the behaviour away from the mode- 
conversion region is given by 

, z oc^ex P (^yX|I d ^ (4.1) 

where (3q is the plasma beta at z = 0. This equation tells us that to obtain a constant 
amplitude in the incident and transmitted waves, as we had in the isothermal case, we 
must make the transformation v z — > v z A 1 / 4 /^^ 2 (shown in the left-hand plot of Figure[5]) . 
Note that if we set A = 1 we revert back to an isothermal atmosphere with the same 
equations as Section [3] 

We may then follow Cairns & Lashmore-Davies to find the behaviour around the 
conversion region, giving exactly the same conversion and transmission coefficients as 
before. So varying the temperature profile has no effect on the amount of transmission 
and conversion as a wave propagates through the mode-conversion region. 

As in the isothermal case, the amount of transmission may be calculated directly from 
the numerical results, so we compare this to the analytical prediction for various values 
of k x . The agreement between the two is again excellent, as illustrated in the right-hand 
plot of Figure 



5. Conclusions 

Using a combination of analytical techniques and numerical computations we have 
studied the phenomenon of mode conversion in a simple, one-dimensional model. More 
specifically, we have studied the downward propagation of linear waves in a gravitation- 
ally stratified atmosphere permeated by a uniform, vertical magnetic field. These waves 
propagate from low- to high-/? plasma, and the mode conversion occurs where the sound 
and Alfven speeds are equal where |3wl. 
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Figure 5. Left: Transformed vertical velocity, v z , for uj = 2-7r\/6 y7/3o/2 and 
k x = vr at t = 7 Alfven times. The dashed line indicates where c s = va- 
Right: Ratio of the slow wave amplitudes in high and low (3, telling us how much of the in- 
cident wave will be transmitted into the high-/? plasma. Here to = 47r\/6\/7/3o/2, the solid line 
is that predicted analytically, and the stars are the values calculated from the numerical data. 

In Section [3] we concentrated on an isothermal atmosphere using a combination of the 
WKB method, and a method developed by Cairns & Lashmore-Davies to examine mode 
conversion both outside and around the conversion region. This allowed conversion and 
transmission coefficients to be found. It was then possible to compare the transmission 
predicted by the analytical method with the numerical data. These showed excellent 
agreement, even as k x becomes large, violating the assumptions. 

This analysis was repeated in Section 2] for a non-isothermal atmosphere. We selected 
a tanh temperature profile to mimic the steep temperature gradient of the transition 
region. Surprisingly the transmission and conversion coefficients found were identical to 
those for the isothermal atmosphere. So the temperature profile has no effect on the 
amount of transmission and conversion as a wave undergoes mode conversion. 

Although the model set-up is not physically realistic, it has given us a great insight into 
the mode conversion process, and the techniques which may be utilised in its analysis. We 
plan to use this knowledge to look at a more realistic model of a two-dimensional coronal 



null point, extending on work done by McLaughlin & Hood (2006) In this situation the 



plasma f3 is infinite at the null point, so as the wave propagates towards the null it will 
indeed be passing from low- to high-/? plasma, as we have examined here. 
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